Isobar model analysis of Tr^ry photoproduction on protons 
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Photoproduction of tt^tj on protons in the energy range from threshold to 1.4 GeV is discussed. 
The data for representative angular distributions recently obtained at MAMI C are analyzed using an 
isobar model. The isobars considered are A(1232) and 511(1535) for Tr^p and rjp states, respectively. 
In accordance with the results of earlier works the main features of the reaction are explained through 
the dominance of the D33 wave with a relatively small admixture of positive parity resonances. 
Comparison with recent experimental results for the photon beam asymmetry is carried out. 
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I. INTRODUCTION 

Several important features of tt'^t] photoproduction on 
protons were explained in earlier papers [IH^] , where dif- 
ferent versions of an isobar model were used as a the- 
oretical basis. According to the experimental results of 
[^-Q the cross section for — n'^'qp exhibits quite a 
rapid rise at low energies and reaches its characteristic 
value a ~ 3.5 at = 1.4 GeV, about 0.5 GeV above 
threshold. Such a behavior might indicate that the cross 
section is mainly governed by the increasing phase space 
so that it is dominated by s-waves. As has been noted 
in [l|, the Born terms seem to be insignificant and the 
reaction is expected to proceed via formation of baryon 
resonances decaying into TrrjN. 

Within an isobar model, production of two mesons may 
conveniently be described in terms of intermediate quasi- 
two-body states, the isobars. These states naturally ap- 
pear if the attractive two-body interaction between the 
final particles lead to narrow resonances. In the case 
of the TrrjN state the ttN and rjN isobars are identified 
with the A(1232) and 5*11(1535) which are well known to 
dominate ttN and rjN scattering. In such a picture the 
observed growth of the total cross section would point 
to the existence of baryon resonances, belonging to the 
third resonance region, which may decay into s-wave ryA 
or TrS'll configurations. One can easily show that only 
in the D33 {P31) state the rjA (TrS'll) system can be pro- 
duced in the relative s-wave. 

The model for Tr'^-qN, with the resonance 

^33(1700) dominating the amplitude, has been devel- 
oped in ref. jj, where ir^r] and tt'^K'^ photoproduction 
were studied within a coupled channel approach. The au- 
thors used the D33 — >■ 77A decay, treated on a tree level, 
as a main driving mechanism. The crucial importance 
of D33 in this calculation is ultimately a consequence 
of quite a large (complex) coupling constant g^A whose 
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value g,,A ~ 1-7 — 1.4j was taken from [7J. The produc- 
tion of the TrS'll state is calculated microscopically as a 
series of interactions 77A — > {-K-qN, ttKA, Tri^TE) — ttSu 
taken up to the first order in the two-body t-matrices for 
rjN and KA (JCE) scattering. 

On the experimental side, recent studies [1, 0: have 
also implied that the resonance I?33(1700) plays a cru- 
cial role in TT^rj photoproduction. In ref. [3| the data 
for 7p — TT^Tjp were explained using the PWA method 
IS] which gave definite quantitative conclusions about 
contributions of different N- and A-like resonances. In 
particular, the authors reported that at higher energies 
{W > 1.8 GeV) the Z?33 wave remains important and is 
presumably populated by the D33(1940) resonance. 

In ref. ^ , in addition to the basic cross section, the 
linear beam asymmetry S was measured. In common 
with the results of [1] the authors of [3] emphasized the 
dominance of the Z?33(1700). The importance of the 77A 
configuration is indicated by the pronounced peak in the 
tt'^P mass spectra around the A mass. The data for the 
E-asymmetry are ingeneral agreement with the predic- 
tions of the model [l|. The I?33(1700) dominance is de- 
duced from the fact that without this resonance one ob- 
tains strong qualitative disagreement with the measured 
asymmetry. 

A remarkable feature of the data reported in ref. 01 
is that the ttSu configuration appears to not be as pro- 
nounced as it may follow from the theoretical calculation 
of ref. [l|. The case in point is that the apparent shift- 
ing of the TT^p peak predicted by the model jH is not 
evident from the experimental results. Instead of this, 
as is shown in Q, the simple model containing only 77A 
decay accounts for the data much better than the so- 
phisticated approach including TrS'll formation. On the 
other hand, trivial assumption that the role of ttSu con- 
figuration in 7p — > rjir'^p is negligible may lead to misin- 
terpretation of the data. It is clear that the S'ii(1535) 
resonance should, in any case, be present on the Dalitz 
plot due to the rjN final state interaction. The reason 
for its 'disappearance' may lie in the fact that the rjN 
scattering cross section does not exhibit a true resonance 
behavior due to the closeness of the resonance mass to the 
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r]N threshold. Furthermore, the presence of an intensive 
■D33 — >■ r]A decay makes it difficult to estimate quantita- 
tively the role of the irSu channel. These observations 
suggest that studying TriV or r]N invariant mass distribu- 
tions does not enable us to clarify the separation of the 
events related to 7r5ii formation. In this connection, it is 
desirable to have methods allowing identification of vrS*!! 
without recourse to the rjN invariant mass distribution. 

The crucial role of the D33 configuration in n^ri photo- 
production is a very important result by itself. It implies 
that, apart from the (7,7r°) and (7,??) reactions deter- 
mined respectively by the A and Sn formation, we have 
another process, [■jjir^r]), whose amplitude in a wide 
energy range is mainly governed by a single dynamical 
mechanism, the excitation of the wave, in particular 
the 1)33(1700) state. Clearly, this feature provides an ef- 
fective method to study the properties of this resonance. 
However, in contrast to single tt" or 77 photoproduction, 
here we face the technical problems associated with three- 
body kinematics, where the particle energies and angles 
arc distributed continuously. It therefore becomes diffi- 
cult to perform a general multipole analysis, primarily 
since there is a multitude of ways to couple angular mo- 
menta of particles to the total angular momentum. 

Important steps on the road to a systematic study 
of the partial wave structure of tt^t] photoproduction 
were made in ref . , where the authors presented a phe- 
nomenological model for — )■ n'^'qN and discussed the 
angular distributions of final particles. Although, in gen- 
eral, unpolarized measurements are not sufficient to dis- 
tinguish between all partial waves, assuming that only 
one of these waves is important, it is possible to reach a 
rather definite conclusion about the corresponding quan- 
tum numbers. The analysis permits the identification of 
the quantum numbers of the dominant resonance states 
simply by analysing the shape of the measured angu- 
lar dependencies, in which the states with different spin- 
parities exhibit their own signatures. It is also empha- 
sized in that some of the observables depend weakly 
on the model parameters. This fact makes such a method 
especially useful. 

The formalism developed in ref. 6] was then applied to 
the analyses of the data recently obtained at MAMI C '3] . 
As a main subject, the authors of [8] studied the angular 
distributions of particles in — >■ Tr'^r/p. This approach 
was shown to be very effective as a tool to study the 
properties of the 1)33(1700) resonance. In particular, the 
ratio of ttSu to rjA decay widths, which, as pointed out 
above, is difficult to extract from the invariant mass dis- 
tribution, was deduced from the data with appreciated 
accuracy. 

Here we continue the analysis of the reaction 7A^ — > 
TT^rjN using the isobar model which was partially de- 
scribed in refs. @ and [1]. The paper is organized as 
follows. In the next section we briefly outline the formal- 
ism, and summarize some results, of ref. @ which pro- 
vides a framework for the calculation to follow. Then in 
Sect. Hill we present our results where the main emphasis 



is put on the angular distributions of the final particles. 
In Sect. nil Bl we consider in some detail the photon asym- 
metry E for 7p — )> n'^rjp and discuss the main properties 
of this observable. In Sec. lIVI we close with a brief sum- 
mary and an outlook. 



II. THEORETICAL BASIS 



We consider the reaction 



(1) 



TT"{uj^,q^) + r]{ujjj,q,^) +Nf{Ef,pf) 

in the overall center-of-mass system. The corresponding 
energies and three-momenta of the particles are given in 
the parentheses. Furthermore, throughout the paper we 
use the following notations for kinematic variables: 



w 

uja, a = ttN, rjN 



photon lab energy, 

total cm. energy, 

invariant energy of the 
two-particle subsystem a, 



= (©i, ^i), i ^ TT,ri solid angle of the momentum 

Qi in the overall cm. frame. 



uj*, q*, J = 7r,?7 



i = TT, ?7 



energy and three-momentum 
of the i-th particle 
in the iN cm. frame, 

), solid angle of the 

momentum q* , related 
to the K{H) system. 



The meaning of K and H systems is explained by Fig. [2] 
(see also text after eq. dS])). 

A theoretical basis of our calculation is presented in 
ref. [1| and partially in ref. We assume that the ttt? 
system does not resonate over the energy region in ques- 
tion. The reaction amplitude is divided into the Born 
and the resonance part. 



(2) 



and is visualized by means of the diagrams in Fig.[T] The 
indices ruf = ±1/2 and A = ±1/2, ±3/2 denote the z- 
projection of the final nucleon spin and the initial state 
helicity, respectively. The Born sector is formed by the 
diagrams (a)-(f). As already noted, its role in the channel 
Tr'^r/p turns out to be insignificant. 

According to the isobar model concept, the resonance 
part of the amplitude may be superposed in the usual 
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FIG. 1: Diagrams representing the amplitude for the 7A'^ -> 
TTT^A'^ reaction in a simple isobar model. The notations A and 
N* are used for the resonances A(1232) and 511(1535). 



fashion, as a sum of two terms ^ 

,R _A-nA) A^N') 



(3) 



corresponding to the diagrams (g) and (h) in Fig.[TJ Each 
term in eq. ([3]) has the form 



R(J^;T) 



xf^l''>!iW,q.,q„Pf), a = rjA,nN*, (4) 

where the summation is over the resonance states deter- 
mined by their spin-parity and isospin T. The latter 
is incorporated through the factor Ct- Since the 77- meson 
has isospin zero, this factor is the same for both the 77A 
and ttN* channels. In our case, where only the states 
with T — 3/2 are assumed to contribute, Ct is 2/3. 

The functions A^{W) are the helicity amplitudes for 
7P R. Their energy dependence was parametrized in 
the simple form 



(5) 



where Mr is the resonance mass and is the cm. pho- 
ton energy, calculated at the resonance position. For the 
dominant resonance, I?33(1700), the exponent ni/2 was 
set equal to 1, whereas the function A^I'^{W) was fitted 
to the data separately in each energy bin in Fig. [5] For 
the other resonances we used n\ =0. 

The factors Gr{W) stand for the resonance propaga- 
tors for which we use the nonrelativistic Breit-Wigner 
form 



Gr{W) 



1 



W-MR + lTg,iW) 



(6) 



^ In the expressions to follow the resonance 5'ii(1535) is denoted 
by AT*. 



In eq. (HI), the functions /^^"\ a = rjA, ttN*, describe 
the decay of the resonance R into TrrjN via intermediate 
formation of 77A and TriV* states. 

We take the Z-axis of the overall cm. frame along 
the photon momentum k. The OXZ plane is spanned 
by the momenta k and Qr, (see Fig. [2]). To fix the final 
state kinematics we choose as independent variables the 
particle momenta associated with the partition rj + (nN) , 
the solid angle of the rj momentum fi^ = (0^,$^ = 0), 
the invariant energy cj^at of nN subsystem, and the solid 
angle of the momentum q* in the ttN rest frame. 

For the data analysis, two types of coordinate sys- 
tems O'x'y'z' , related to the ttN cm. frame, are used. 
In the first one, the canonical (K) system (Fig.[2][a)), 
the axes are chosen parallel to those in the main cm. 
frame OXYZ fixed in space. In the helicity (H) system 
(Fig. [21b)) the z'-axis is taken to be along (f^ +pf = —q^. 
For a single event the systems are connected via appro- 
priate rotation by the angle tt — 0^ around the j/'-axis. 
After integrating over G,, is carried out this equivalence 
does not hold, and in this sense the distributions W{nf) 
and W(n^), i = IT,?] are independent. The values re- 
lated to the K and H frames are further denoted by the 
indices K and H respectively. 

In the K system, the pion momentum, q*, is related 
to the pion momentum (f^r in the overall cm. frame via 
the Lorentz boost along the vector —9,,, giving 



qn ^ q* + Xj^qrj , Xj, 



(^7tN + W -UJri 



(7) 



In Sect. lIIIBj we also use the partition tt -|- {riN), tak- 
ing as independent kinematical variables the solid an- 
gle ri^ = (Ottj^tt = 0), the invariant energy uJrjN, and 
the spherical angles of the rj momentum in the rjN cm. 
frame. In this case the reaction plane is spanned by the 
vectors k and q^. The relation between the momentum 
of the 77 meson is similar to that in 



q-n = q., 



n 



Xr^qT, 



Xri — 



(8) 



We adhere to the non-relativistic concept of angular 
momentum and do not include, for instance, effects of 
spin transformation caused by the Lorentz boost. Then, 
in the coordinate system chosen, the angular dependence 
of the amplitudes may be decomposed by means of 
spherical functions. In the K system, associated with the 
partition rj + {ttN) (Fig. [2]) we will have 

rR{riA)K _ r^R{r,A)(,, nV^^|Ma 

X ^I'm. fM. Yu^^) 4;,o(0.) . (9) 

X E Ct-ZM,-™ ImYU^^) 4/:-™o(0.) • (10) 
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(a) 



FIG. 2: The coordinate systems {x'y'z') used for the analysis 
of angular distributions of pions in the np rest frame. In the 
canonical system {K system) the z' axis is taken parallel to 
the beam direction, whereas in the helicity system {H system) 
it is aligned along the total irp momentum. In both systems, 
the x' axis is in the reaction plane and the y' axis is chosen 
as y = (p^ X kj)/\p^ x kj\. 



cients for the coupling ji 
determined as 



are the usual Clebsch-Gordan coeffi- 
j2 = js . The coefScients Ai , 



TABLE I: Angular momenta associated with a decay of the 
resonance R(J^) into the channels 7)A (denoted by L,,) and 



J^{L2T2.]) 






r(^3i) 


2 


1 


r(j^3i) 


1 







0,2 


1 


r(j^33) 


1,3 


2 


§'(1^35) 


2,4 


3 


§'■(^35) 


1, 3 


2 



tions normalized to unity. For instance, 



_ (X^q,X / {2L^-l){2L^)\ 
' \ ql ) \j {2l-l){2L,-2l)\{2l)\' ^ ' 

stem from the expansion of the function Yl^a/^ {^-n) over 
the products of the spherical functions depending on il^ 
and r2^. 

The corresponding expressions for /^^"''^ can easily be 

obtained from eqs. © and (ITUl) via rotation of Yim{^^) 
by the angle 0,, around the Y axis. 

The coefficients a = r/A, ttA^* are 

parametrized in terms of coupling constants 

f«(-^*)K^) = /^-^-{^''^^ G^.K^)g;^', (13) 

where Ga and Gat. are the A and N* isobar propagators. 

Unlike ref. [l[ here we use a phenomenological rather 
than a dynamical approach to describe the transition to 
■kN* . Namely, the vertex R — >■ ttN* appearing in the 
diagram (h) in Fig.[T]was parametrized in the form of an 
ordinary resonance decay R — ^ 'kN* fixed by the constant 
/rttN' (see eq. (fT^ ). In contrast, in ref. [l| only the 77 A 
decay is treated as a 'primary process' whereas ttN* is 
generated during the rjp final state interaction. As noted 
in 0, we expect our approximation to be reasonable at 
least around the resonance maximum. 

The angular distributions are analysed in terms of the 
functions W{6) and W{4>) which are differential cross sec- 





W{<t>^)^l /^sine^e, (15) 

^ 

with a standing for the total cross section. 

To effectively include the total cross section into the 
fitting procedure we use the products d-W{0^^^'') and 

aW{4>^^^''), where a is the total cross section, averaged 
over the corresponding energy bin. 

The adjustable parameters are collected in Tablellll 
Since the resonance amplitudes are proportional to the 
product of the electromagnetic and hadronic couplings, 

the amplitudes Ax and the ratios /?„ = ri-'^Ar/^'tot' ^ — 
rjA, ttN*, cannot be well determined individually. For 
this reason, we use the quantity \/A7A^i/2 together with 
the ratio of the helicity amplitudes and A1/2, and 
the ratio r of i? — >■ ttN* and R — >■ rjA branchings. 

We note that our fitting procedure did not include 
varying the total widths. The reason lies in the close- 
ness of the resonances, especially i?33(1700), to the ttij 
production threshold. As a result, the set of angular 
distributions alone does not impose sufficiently stringent 
constraints on the values of the resonance widths, and 
to some extent, on their masses. In this respect, we fix 
the widths, taking their values from the PDG compila- 
tion [l^] or from the references cited there (the width of 
-D33(1700) was taken from the recent analysis of [ll| and 
that of P3i(1750) from ref. [HI). For the same reason, 
the masses of the well known resonances, rated by four 
or three stars, were varied around their values given by 
PDGtlOl]. 
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FIG. 3: Powers of 77A and nN* relative momenta averaged 
over the isobar distribution as functions of the total cm. en- 
ergy. As indicated in the legend, the curves determine energy 
dependence of partial decay widths of the corresponding reso- 
nances (see eqs. (|16|l and HZ}). As noted in text, for R — rjA 
decay only the lowest value of the orbital momenta Lr, is con- 
sidered. 



The orbital angular momenta and related to 
the r]A and ttN* decays of the resonances are listed in 
Table m As discussed in Sect.HIEl in the D33(1700) ^■ 
77A decay we included both s- and d-waves. As for other 
resonance states, for the sake of simplicity, we take only 
the lowest of the two possible values of L,^ (see Table|l}. 



III. DISCUSSION OF THE RESULTS 

The general properties of the measured angular dis- 
tributions and energy spectra were considered in ref. [1] . 
To interpret the data, the authors of ref. P used the sin- 
gle resonance model. It was assumed that there is one 
resonance (or, in general case, one partial wave) domi- 
nating the amplitude in a wide energy region. Following 
this concept they retain only one term in the sum ([4]) 
to roughly reproduce the main features of the observed 
angular distributions and thus to reach a tentative con- 
clusion about the quantum numbers of this dominant res- 
onance. 

The results obtained in ref. Q are now used as a start- 
ing point for a more refined analysis. Namely, having 
established the appropriate magnitude of the dominant 
partial wave, we successively add other waves in order to 
reproduce the structural details of the measured quan- 
tities. Thereupon we calculate the angular distributions 
measured in ref. as well as the photon beam asym- 
metry presented in Q and fl3| , and compare our results 
with the data. 



A. Angular distributions 

Before dealing with the results of our fit, it is instruc- 
tive to qualitatively discuss the role of different partial 
waves in the angular distributions obtained in the ex- 
periment. Firstly, we consider the near-threshold region 



where the energy is not high enough to produce iso- 
bars in higher angular momentum states, so that spin- 
parity selection rules may be used at least for qualita- 
tive estimations. In essence, the question is related to 
the energy dependence of the partial decay widths for 
i? — > a — !■ TrrjN, a — 77A, ttA^*, which in turn are de- 
termined by the appropriate powers of the rjA and TriV* 
relative momenta averaged over the isobar distribution 
in the available phase space: 



2L^ + 1 



Lo duj . 



(w-MA)2+rA(c^)V4 

W-m 

.2L. + l(W/,c^) 

Mjv +mn 

rAr-(a;) 



(16) 



vv — m^r 



(L^-MAr-)2 + riv*(a;)2/4 



uj duj . (17) 



The particle momenta appearing in eq. (jl6p and (|17p 
are determined by the usual kinematic formula 



q^{W,u)^ ^ ' * 



n,r] 



(18) 



with A standing for the triangle function A(x, y, z) — {x — 
y - z)"^ - iyz. 

In Fig. [3] we present the ratio 
(g|^^'+^)(W^)/(gf^*+^)(M), i = 7r,r/, calculated for 
M = 1.72 GeV with typical isobar masses and widths: 

Ma = 1232 MeV , Fa (Ma) = 120 MeV , 



1535 MeV , Fa,. (Mat. ) = 150 MeV , 



As might be expected, in each individual channel, 77A or 
ttN* , the energy dependence of the partial cross sections 
exhibits the typical threshold behavior ^ 
i = -K, rj. This simple relation, however, is slightly vio- 
lated in the real situation when both configurations, r/A 
and ttN* , are included. Firstly, because of the different 
parities of A and N* the values L,, and L^, entering the 
same partial wave J^, always differ by 1. Secondly, the 
nominal thresholds corresponding to the production of 
stable A and N* are = 1.22 and 1.01 GeV, respec- 
tively, for rjlS. and TriV* channels. As a result, the r/A 
channel is effectively opened at energies about 200 MeV 
higher than those associated with ttN* production. Fur- 
thermore, the smaller mass of the pion, in comparison to 
the ry, leads to the availability of a larger phase space for 
the ttA'^* configuration already at rather low energies. By 
referring to Fig.|3]it may be seen that in the 77 A channel 
the I?33 state exhibits the largest rate of growth. The 
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FIG. 4: Estimate of resonance contributions to "/p — >■ Tr^r^p 
according to formula ([3]). Presented are the resonances which 
were examined in the quark model analysis of ref. [ij as the 
baryons decaying into the rjA channel. The dashed lines show 
the contribution of P33(1600) and P33(1920). For the electro- 
magnetic and hadronic decay widths the values from the PDG 
compilation ilO] are taken. For the partial widths T-nrjN the 
same relation r,rj7jv = 0.1 Ttot for each resonance is used. 



next most important waves are P31 and P33, producing 
the 77 A system in the relative = 1 state. Further- 
more, the P31, leading to the s-wave in the ttN* mode, 
may kinematically be the most favorable wave close to 
the threshold. 

Some A resonances, which may be coupled with the 77 A 
state, were calculated within the constituent quark model 
of ref. [1^ . A rough estimation of their contribution can 
be made by taking the corresponding parameters from 
the Particle Data Tables [l^ . Neglecting background and 
using unitarity one arrives at the sum over the resonance 
states R 



with the partial cross section reading 



TT 



a^{W) = —{2jR 



1) 



{w-MRr + rgfiw)/4 



(19) 



(20) 



The cross section predicted by formula is plotted in 
Fig.|4] together with the contributions of the individual 
resonances. For the partial decay widths we used the 
common relation TT^riNiM) — O.lTtotiM) for all states. 
Although this estimate is very rough because of the un- 
certainties in the partial widths Tt^jjN, one can expect 
that the P33(1600) resonance already mentioned can re- 
ally be visible near the threshold, whereas the P3i(1750) 
and P35(1905) may come into play at higher energies. 

Since we expect the D33 to dominate, it is useful to 
have an analytic expression for the cross section contain- 
ing only -D33. Retaining in the sum Q the term with 



J-^ = 3/2", and using in eqs. ^ and ^ L,, = (for 
the moment we omit the d-wave in the D^^ — > t^A decay) 
and = 1 from Table U we obtain in the H system 

da 



duij^Ndfl^ d cos 

^ (a?/2 + Al/^) (|ci|2 + \c2f + 2Re (cics) cos^f 



2 l^V2 



A 



3/2 



- (1-3 cos" e„) 



X (|cip + |c2p + 2Pe(ciC2)cos0f - -jcipsin^^^ 
+ (^2|cip cos 6if -(- Pe(ciC2)) sin2e^sin6lf cos(?!)f 



Icipsin^ sin^ 9^ cos2(/)f 



(21) 



where we omitted unessential kinematic factors. The co- 
efficients ci and C2 in ((2T|) depend only on cos 9^ and 
cJttTV, and are independent of 0^ and ip^ . In terms of the 
functions P^33(a) ^ and ^ they read 

Ci = (^i^i533('7A) _^^Z)33(^A^*)^ ^* ^ (22) 

C2 = P^-("^')x,g^. (23) 

The factor in eq. comes from the Lorentz trans- 
formation to the ttN rest frame q*^ — >■ q* (see eq. ([7])). 
The angular distribution W{^1^ ) can easily be obtained 
from (|2ip via appropriate integration. Since ci and C2 are 
independent of 0,,, the integration over cos Qri is trivial 
and we arrive at 



Wi9^)^ A + Bcos9^, 



(24) 



with 



A = — 
N 



(|cir + |c2r)dw^Ar 

B = ^ I Re{ciC2)dujT,M , 
where the normalization factor reads 



2 

N 



(25) 



N^J (|ciP + |c2p 

+ 2Re{ciC2)cos9^)duj^Ndcos9^ . (26) 

It is worth mentioning that the quantities A and B in 
eq. are still functions of cos9^ . 

The distribution over the azimuthal angle 0^, ob- 
tained via integration of the cross section (|2ip over 
cosO^, cos9^ , and w^tvj has the form 



WicP^) = ^(l+^—^Ccos2<p^ 



(27) 



The coefficient C is independent of cj)^ and reads 

C^^y \ci\HmH^ duj^Nd cos 9^ . (28) 
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FIG. 5: Angular distributions of pions calculated in the canonical and helicity coordinate systems (the first four columns) and 
angular distribution of rj mesons in the overall jp cm. frame (the fifth column). The data are taken from @]. The curves 
are our isobar model calculation corresponding to the solutions presented in Tablellll solution I (solid), solution II (dashed), 
solution III (dash-dotted). The dotted curve is obtained within the single model, including both D33(1700) and D33(1940) 
resonances. In the last case the parameters of both resonances were taken from the solution I. 



The parameter a is determined by the ratio of the D33 
helicity amplitudes as 

a{W) = . (29) 

Several conclusions important for understanding the 
experimental data presented below can be drawn. Firstly, 
as is found in ref. |8j, the distribution W{9^) is rather 
sensitive to the relative contribution of the ttN* channel 



to the 1)33(1700) decay. Quantitatively this sensitivity 
was discussed in |8| in terms of the parameter 



-D33— >-7rjj7V 

where T^j^^^^^, a = rjA, ttN* , stands for the corre- 
sponding partial decay width. 

From ([23|) and (f24|) one sees that if the contribution of 
ttN* vanishes, W{9^) becomes isotropic {B = 0). The 




FIG. 6: Ratio ^3/2/^1/2 of the helicity amplitudes for 
£'33(1700) as a function of the lab photon energy. The solid, 
dashed, and dotted curves correspond to the solutions I, II, 
and III from Tablellll The empty circle shows the ratio ob- 
tained in ref. [ill], whereas the filled circle is the value advo- 
cated by PDG 



addition of the nonzero — TriV* coupling results in an 
approximately linear dependence on cos 0^ (there is an 
additional dependence of the factors F^^^^^-'^ ^ and 
in (22) on cos 9!^, which is, however, rather weak). The 
corresponding coefficient B is proportional to Re{ciC2) 
p4)) and is therefore rather sensitive even to a small ad- 
mixture of TriV*. 

In the canonical frame one has 

da 



dujTrNd^n cos 8,, 

^?/2 ( I ci P + I C2 P - 2i?e (ci C2 ) cos 61^^ ^ 

+ (^3/2 - ^?/2) (iciP sin^ 0f + |c2p sin^ 6^ 
+ 2i?e(ciC2) sin 9^ sin 9^ cos 0f ) , (31) 



with 9Trrj being the angle between q* and g^. Although 
the formula ([3T|) looks less cumbersome than eq. (|2T]) , the 
coefficients ci and C2 are now functions of all four vari- 
ables uJt^n, J ^^'^ cos0^. Therefore, a simple general 
analysis, like in the case of the H system, is not possible. 
The only conclusion we can draw is that since ci and C2 
depend rather weakly on the angle cj)^ , the distribution 
W{(l)^) can approximately be represented in the form 



1 

2^ 



1 



(32) 



where the coefficient C is equal to 



C = - / i?e(ciC2) 



(33) 



X sin 8^ sin 9^ du^N d cos 6r, d cos 9^ 



A look at eqs. (p3)) and ([27]) shows that, contrary to 
W{(j)^) whose form depends only on the parameter a, 
in the K frame the (?i^-dependence is solely due to inter- 
ference between 77A and ttN* mechanisms. We therefore 
expect that the shape of W{(j)^) is much more sensitive 
to the ttN* contribution than the 07r-distribution in the 
helicity system. This assumption is well confirmed by the 
analysis of ref. Q . 

In Fig. [5] we demonstrate our results for different an- 
gular distributions. As noted above, we choose as in- 
dependent kinematic variables those associated with the 
Tj + (nN) partition. Namely, the results are presented 

for the distributions W{9^'^"^) and IF(0f ^^^) over the 
pion angle in the K and H system (see Fig. [5]), as well 
as for the distribution M^(8^) over the angle of the rj 
momentum in the overall cm. frame. The solid, long- 
dashed, and dotted curves present three different solu- 
tions, providing an acceptable description of the data. 
The resulting parameter sets corresponding to each solu- 
tion are collected in Tablellll The partial cross sections 
in the last column in this Table are calculated according 
to formula (l20l) at the energy W = AI^, and give ap- 
proximate contributions of the corresponding resonances 
to the total cross section. 

It is worth mentioning that the fits collected in Tablellll 
do not necessarily give the lowest value of x^- Rather, 
we present them as three sample results which describe 
the same data via different dynamical contents. Namely, 
solution I (first row in Tablellll solid line in Fig. [5]) is 
mostly consistent with our assumption about -D33 dom- 
inance over the whole energy region considered. Here 
other partial waves are of relatively minor importance. 
Solution II (second row, long-dashed line in Fig.[S|) de- 
scribes the data at > 1.2 GeV via an appreciably 
large contribution of the resonance ^33(1920). Finally, 
solution III (third row, dotted curve in Fig.[S|) contains a 
small admixture of the d-wave in £>33(1700) -> 77 A decay, 
which is neglected in the first two cases. 

The following conclusions can be drawn: 

(i) The general properties of the data in Fig.[S] (the 
number of minima and the convexity direction) are 
roughly described with only the D33 wave. This fact, 
considered as a strong evidence of the I?33 dominance 
was already discussed in ref. 0. In particular, one can 
see that at energies above = 1.2 GeV the experi- 
mentally observed shape of W{9^) and W{(f)^) is rather 
close to the predicted form a + b cos 9^ and c + d cos 2ip^ 
(eqs. (IM|) and (|27p '). In the canonical frame the distribu- 
tion W{(p^) is also in good agreement with the analytic 
form A + B cos (j)^ provided by the single D33 model (see 
eq. (|32p ). As is shown in ref. Q, other partial waves fail 
to describe all the angular distributions simultaneously. 

(ii) Since the coefficient C in eq. (I?7|) is a constant, the 
(j)^ dependence is totally determined by cos 20^. Then, 
the convexity direction of W{(j>^) (up or down) is gov- 
erned by the sign of the difference 1 — a. Therefore, the 
observed cj)^ -dependence allows one to make a conclusion 
about the value of the parameter a. More specifically, ac- 
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TABLE II: Parameters of the isobar model fitted to the data in Fig. [5] The first three rows give the parameter sets, to which we 
refer in the text as solutions I, II, and III. The forth row lists the corresponding numbers given by the Particle Data Group 
The sign of ^ is the sign of the product of A1/2 and the R — >■ r/A coupling constant /aatj (see eq. (HU). The relative 
sign between R —>■ rjA and R — >■ ttN* couplings is taken positive for all resonances. The partial cross sections are calculated 
according to the approximate formula (|20|) at = Mr. An extremely large (small) value of the ratio r (see, e.g., those for 

P3i(1750) and ^35(1905)) means that the corresponding partial decay width F^''^'^' ') is close to zero. X^/^ = 2.76, 3.06, 

and 3.39 for the solutions I, II, and III, respectively. 



J"[L2T2j(Mfl)] 


Mr 

[MeV] 


rf.t(A/B) 

[MeV] 


^/^A^f2{MR) 
[10-3GeV-i/2] 






a'' (Mr 
[Mb] 




§-[1)33(1700)] 


1701 ± 1 


375 


10.6 ± 0.2 


0.95 ±0.01 


0.70 ±0.03 


0.186 






1720 


375 


15.1 


0.95 


0.67 


0.384 






1709 


375 


14.8 


1.00 


0.67 


0.370 






1700 


300 




0.82 








f+[P33(1600)] 


1657 ± 5 


350 


-0.10 ±0.011 


0.39 ±0.06 


1.78 ±0.32 


6.70 ■ 10" 


6 




1636 


350 


0.17 


1.20 


169 


7.02 ■ 10" 


3 




1670 


350 


0.18 


3.71 


71.4 


1.90 • 10" 


2 




1600 


350 




0.39 








i+[P3i(1750)] 


1832 ± 1 


400 


8.35 ±0.9 




0.22 ±0.11 


3.78 • 10" 


2 




1820 


400 


9.93 




2.6 ■ 10"^ 


4.53 ■ 10" 


2 




1826 


400 


40.2 




1.0 ■ 10"^ 


0.720 






1750 














|+[P35(1905)] 


1873 ± 4 


330 


-25.5 ±0.6 


-0.70 ±0.03 


(9.1 ±7.5) ■ 10"^ 


0.506 






1890 


330 


-10.8 


-1.64 


1.7 ■ 10"^ 


0.231 






1890 


330 


9.3 ■ 10"" 


-1.94 


1.1 ■ 10^ 


2.38 ■ 10" 


3 




1890 


330 




-1.73 








|+[P33(1920)] 


1894 ± 3 


200 


11.87 ±0.42 


1.15 ±0.06 


0.14 ±0.09 


0.319 






1842 


200 


15.5 


1.49 


0.89 


1.29 






1848 


200 


4.42 


1.62 


6.4 


0.462 






1920 


200 




0.575 








|-[P>33(1940)] 


1870 ± 1 


450 


19.9 ± 1.1 


1.65 ±0.02 


0.22 ±0.03 


0.695 






2069 


450 


-37.6 


1.77 


0.10 


2.26 






1937 


450 


-7.8 ■ 10"^ 


5.09 


1.1 ■ lO*' 


1.04 






1940 
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T3 
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COS 0, COS 0„ COS 0„ 

FIG. 7: Our prediction of angular distributions for energy 
bin 1.7 < W < 1.9 GeV compared with the CB-ELSA data 
of ref. [3] (open circles). In the first row, the angle 9ij is the 
angle of the ith particle in the ij cm. frame with respect to the 
momentum of the third particle. In the second row, O is the 
angle of the corresponding particle in the overall cm. frame. 
The solid (dashed) curve is obtained with the parameter set 
I (II) from Table HIl The dotted curve includes only the D33 
wave. The theoretical results have arbitrary normalization. 



cording to the expression (P7|) . the characteristic shape 
of W{(j)^) with the minima at <f)^ = 7r/2 and 0^ = 37r/2, 
visible at > 1.2 GeV, should point to the fact that at 
this energy a > 1 (eq. ([27])). This result is also confirmed 
by W{6^) (first row in Fig. [5]) which, as follows from ([5^ 
and (j33p , would otherwise be convex down in contrast to 
the data. 

Our results for the ratio ^3/2/^1/21 a function of 
the photon energy, are presented in Fig. [51 As may be 
seen, version III requires an unexpectedly sharp energy 
dependence of this parameter, so that already at Ej = 
1.4 GeV it amounts up to 2.4, thus visibly exceeding the 
average PDG value 0.817 ± 0.242. At the same time, 
versions I and II allow one to suppress the undesirable 
increase of this ratio at higher energies. 

(iii) According to our qualitative discussion above (see 
text after eq. ([501) ). the visible forward-backward asym- 
metry in W{6^) may be interpreted as a signal of a 
nonzero ttN* contribution. This conclusion is trivial in 
the sense that the ttN* channel always contributes to 
the reaction — > ir'^rjp via rjp interaction in the final 
state. However as noted in the Introduction, direct mea- 
surement of the rjp spectrum to study the role of this 
channel may fail, primarily because of the overwhelming 
contribution of the 77A decay mode. At the same time, 
the role of TriV* may be estimated rather precisely using 
the asymmetry rate of W{6^). Indeed, as is shown in 
0, here the different states J'^ are added incoherently, 
so that the shape of W{6^) is always symmetric if only 
the channel rjA is taken into account. In this respect, 




E, [GeV] [GeV] 



FIG. 8: Total cross sections for jp — )■ rj-K p. The the- 
oretical results are obtained with parameter set I in Ta- 
ble [III Left panel: full calculation (solid curve), -D33 wave 
(dashed curve), P33 wave (dotted curve). Born terms, shown 
by the diagrams (a)-(f) in Fig.[Tl (dash-dotted curve). Right 
panel: the contribution of 1)33 (1700) and D33(1940) (dashed 
curves), ^35(1905) (solid), P33(1920) (dotted), P3i(1750) 
(dash-dotted). 



the observed asymmetry is a direct consequence of the 
interference between 77 A and vrA^* mechanisms. 

(iv) In the region E^ — 1.2 — 1.4 GeV, the distribu- 
tion W{&rj) over the rj angle in the overall cm. system 
demonstrates the characteristic shape a + 5cos^ G, van- 
ishing at lower energies. This feature is explained in dif- 
ferent ways in our fits I-III. For example, the parameter 
sets I and 11 describe 11^(8,,) via a relatively strong con- 
tribution of ^35(1905) and P33(1920) respectively . In 
the third version of the model, the increasing role of the 
cos^ Q term is ascribed to the d-wave fraction in the r/A 
decay of 1)33(1700). Since the D33 mass is pretty close to 
the TTT] threshold, even a small d-wave admixture in the 
77A mode may lead to a sizable contribution of this wave 
at energies above E^ — 1.2 MeV. The dotted curve in 
Fig. [5] is obtained by assuming that the d-wave fraction 
comprises only 2% of the total D33 — 7r7]N width. 

(v) In Fig. [7] we show our prediction of some angular 
distributions for three different partitions, tt" -|- (rjp), ri + 
(tt'^p), and p + {tt'^i]), together with the data obtained in 
ref. 0]. As one can see, our calculations are in general 
agreement with the experimental results, perhaps except 
for the distribution over cosfl^p. At the same time, these 
observables demonstrate a rather weak sensitivity to the 
model ingredients. Already the single model (short- 
dashed line in Fig.[71) correctly predicts the main features 
of all distributions. Addition of other resonances does not 
visibly change the calculation. On the whole, all three 
parameter sets in Tablellll give equally good descriptions 
of the data. 

(vi) As noted in the Introduction, the wave seems 
to also remain important at higher energies. We base 
this conclusion on the fact that the single D33 model 
works quite well at E^ > 1.3 GeV (short-dashed line in 
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Fig. [5]). At the same time, our prediction for the total 
cross section with only Z?33(1700) starts to already un- 
derestimate the data at = 1.1 GeV. This may be 
considered as an indication of an additional source of the 
D33 wave which should come into play at higher energies. 
In ref. [3| this effect is ascribed to 1)33(1940) which is clas- 
sified as a one-star resonance in the PDG listing fld\. In 
the present work, to maintain the growth of the partial 
cross section a^^^ we also introduce this resonance with 
free parameters. However, one has to keep in mind that 
since 1)33(1940) is at the upper limit of our energy region, 
the fitting results are not very sensitive to the exact con- 
tribution of this resonance. The same must hold true for 
^33(1920). 

(vii) Fig. [5] demonstrates our calculation of the total 
cross section, where contributions of individual terms en- 
tering the amplitude are plotted separately. The results 
are obtained with the parameter set corresponding to the 
solution I in Table|lll We prefer this version of our model 
since, as will be shown in the next section, it also yields 
a reasonable description of the beam asymmetry, in con- 
trast to the solutions II and III. The reader should not be 
surprised by the quality of our description in Fig. [S] since, 
as noted in Sect.lIIl the total cross section was effectively 
included into the fit (see the comment after eq. (|T5|) ). 

As predicted by solution I, the major fraction of the 
reaction yield is provided by the coherent sum of the res- 
onances /?33(1700) and i?33(1940) (short-dashed line on 
the left panel in Fig.©. The contribution of 1)33 (1940) 
is comparable to that of 1)33(1700). As pointed out 
above, the need for the former resonance in our calcu- 
lation is mainly dictated by the necessity to prevent the 
1)33 wave from dying out when 1)33(1700) rapidly falls-off 
above > 1.2 GeV. Otherwise the model leaves room 
for P33(1920) leading to solution II. As follows from Ta- 
bleini in this case the contribution of P33(1920) amounts 
up to 40% of the total cross section at Ej = 1.4 GeV. 

It is quite remarkable that the role of ^33(1600) in 
our analysis turns out to be insignificant. The reason 
is that in the near-threshold region some fraction of the 
p-wave comes from the background. The corresponding 
p-wave admixture turns out to be sufficient to obtain a 
more or less accurate description of the forward-backward 
asymmetry in the angular distribution W{Q,-i) of 77 me- 
sons in the overall cm. system. At the same time, at 
higher energies, where the role of D33 strongly increases, 
the background is unable to give the necessary p-wave 
fraction. As a result, our fit requires inclusion of other 
resonances with positive parity, in particular P3i(1750), 
F35(1905), as well as P33(1920). 



B. Beam asymmetry 
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FIG. 9: Beam asymmetry for 7p — >■ Tv^rip averaged over three 
intervals of the beam energy. The data are from ref. [J. The 
solid (long-dashed) curve is our prediction corresponding to 
the parameter set I (II) in Table HTl The short-dashed curve 
includes only the D33 wave. 



this observable demonstrates several properties, which 
we would like to discuss. Firstly, the data show (see 
Figs.ini and [TU|) that in the ir'^p case, when the reaction 
plane is fixed by the condition 0,, = in the overall cm. 
system, the value of E (hereafter denoted by I](7riV)) is 
comparable to zero in almost the whole region of uj^^^, 
except for the lower boundary of the spectrum. In the 
channel ryp, where (p^r — 0, the beam asymmetry 'E,{riN) 
is negative and its absolute value comprises about 0.5. 



The beam asymmetry S for linearly polarized photons 
has been measured in refs. Q and [l3|. It was presented 
as a function of the invariant mass of three two-particle 
subsystems n'^p, r/p, and Tr'^ry. In the region Ery < 1.5 GeV 



Both features may be qualitatively explained by the 
dominance of the 1)33 partial wave. Indeed, keeping only 
1)33 in the sum Q , for the polarized cross section one ob- 
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FIG. 10: Same as in Fig.|9]but for other energy bins. The data 
are taken from ref. 13] . The sohd (dotted) curve corresponds 
to the solution I (III) in Table |IT] 



tains (omitting unessential angular independent factors) 



, da 

dui -jrr, 



1 f da 

2 \da7 



da\ 



doj 

/2^3/2 



^1/2^ 



V3 



sinSf dd^ 







(|cipi(3cos^e-l) + |c2p 
25Re(ciC2) COS 6if) , (34) 



where the coefficients Ci are determined in eqs. ((22)) and 
(|23l) . In eg. (|34l) a± {a\\) is the cross-section with the 
photon beam polarized perpendicular (parallel) to the 
reaction plane. The corresponding unpolarized cross sec- 
tion reads 



da if da\ 



daw 



duJT, 



2 Y duj-jYp diu-^p 



2(^/2 



2 , ^2 



3/2 



sm 



X (|ci|2-f |c2p + 23?e(ciC2)cos0f) . (35) 

It is now straightforward to obtain an analytic expression 
for the beam asymmetry dividing eq. (|34|l by eq. (|35|) . 



In the partition tt -|- (rjN) the expressions for S ^^"^^ 



and for the unpolarized cross section 



da 

duJriN 



are formally 



identical to that given by (l34l) and ([35 
„ da 



doj. 



rjp 



V3 



sin6l,f d6i^ 





2 aH 



,3cos^ 9' 
2^e(bib2)coae^ 



1) 



(36) 



dcr 1 



diij„ 



^(^V2 + ^3/2) / sm<< 



X ( I61P + I62P + 23*6(6162) cos^^j , (37) 



except that the integration variable is now and the 
coefficients q are replaced by bi, determined as 

bi = (38) 

62 - (^F^ss(riA)^^ _ _^ ^D33(^W*)^ (39) 

With increasing cja, the value of Y,{7tN) rapidly goes 
to zero, whereas 'E,{riN) demonstrates a rather smooth 
dependence. This qualitative difference may be explained 
by the kinematic effect caused by the rather large mass 
difference of the tt and 77 mesons. 

Indeed, because of the small mass of the pion, its mo- 
mentum in the major part of the reaction phase space 
is essentially higher than the 77 momentum. As a con- 
sequence, in the partition 77 -I- (ttN) we will have (see 
eqs. (EH) and (^5)) ) 



2 / * \ 2 

»1. 



(40) 



For example, taking for the parameter r the value r = 
0.7, as predicted by the fit I, we will have |ci/c2p « 27 at 
W = 1.7 GeV, which further increases to 74 at W — 1.9 
GeV. 

Since in the cross sections (|34|) and ([35]) the coefficients 
Ci appear in the integrand, it is more appropriate to com- 
pare directly the corresponding integrated values. As an 
example, in Fig.[TT]we show the integrals 

7r 

= y Icpsin^f d0f , i = l,2, (41) 


J, = J \b,\Hme^ d9^ , * = 1,2, (42) 


calculated at Ej = 1.2 GeV, as functions of the energy 
uja, ct = rjA, nN*. As noted above, in the region of very 
low WjrTV where g* ~ we have Ii <^ l2- However, with 
increasing WjrJV the pion momentum q* rapidly increases, 
so that the relation Ii ^ I2 holds in almost the whole 
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FIG. 11: Integrals l|4ip and (|42|l calculated at photon energy 
£■.^ = 1.2 GeV. The values are given in arbitrary units. 



region of lj-hn- For instance, as may be seen from Fig. 111! 
already at w^at = 1.15 GeV we have /1//2 — 10^. 

At the same time, in the polarized cross section 
the coefficient |cip enters with the weight P2{cosd^) — 
(3 cos^ 6^ — l)/2 (see eq. which strongly reduces its 

contribution after integration over 0^ . This last effect is 
due to the rather smooth dependence of |cip on cos9^ . 
As a consequence, the value of I]{ttN) turns out to be 
roughly equal to 



E(7riV) 



(43) 



In other words, in the region where the I?33 wave domi- 
nates, the nonzero value of the asymmetry T,(ttN) is due 
to the nontrivial contribution of the ttA'^* configuration 
to the I?33 decay. For this reason, the rapid decrease of 
T,{ttN) at moderate cOttn is due to the decrease of this 
contribution to the integral in eq. (|34|) . The latter has a 
kinematic reason and, as we have just explained, is even- 
tually caused by the mass difference between tt and 77 
mesons. 

In the second case, when the reaction plane is fixed 
by 0^ = (partition rj + (nN)), the situation is just the 
opposite. Namely, (see eqs. (|38 | - (|39)) ) 



so that 



\^ivN)\ 



^1 «1, 



J2 



J1 + J2 



(44) 



(45) 



In the same limiting case uJrjN — ^ (61 — )■ 0), we get 
from cqs. and (1571) a simple relation between the pa- 
rameter a ([29]) and the value of T,{'riN) 



V3 1 



(46) 



Taking a 
region 



1.2 (the average of this parameter in the 
1.0- 1.4 GeV) we win have I](?7iV) = -0.43, 



which is consistent with the experimental result. Since 
|S(77iV)| dUl) has maximum \^{r]N)\ = 1/v^ at a = 1, 
it appears to be rather insensitive to a as long as a does 
not strongly differ from unity. 

In a similar manner, at lottn — >■ we have ci ~ g* — ?■ 
and the value of Y,{ttN) depends only on the parameter 
a (j29l) . As a result, in the single Z?33 model (or generally 
in the single resonance model) one has the simple relation 



E(7riV) 



(47) 



It is also worth noting that, according to the formula 
([M]). the negative value of the asymmetry Y,{riN) is a 
signature of the fact that A^^^^ and ^3/2^ should have 
the same sign. 

Our results for the photon asymmetry predicted by the 
three solutions discussed in the previous section are pre- 
sented in Figs. [Ul and [TUl As we can see, only version I of 
our model where the 1^33 wave dominates in the whole 
energy region is able to reproduce the main features of 
the measured asymmetry. As would be expected, the pre- 
dicted behavior is close to that given by the single 
model, so that inclusion of other resonances does not give 
rise to a crucial change. At the same time, version II con- 
taining an appreciable contribution of ^33(1920) (see the 
last column in Table llH) yields a rather poor description of 
Y,{riN) in the region 1.3 < i?..^ < 1.4 GeV. As one can see 
from Fig.[9l at low invariant energies LUrjN this resonance 
causes an undesirable shift to positive values, in contrast 
to the data demonstrating a rather smooth dependence 

on UJrjN- 



IV. CONCLUSION 

We performed a phenomenological analysis of the 
photon-induced Tr'^ry production at energies within the 
region 1.0 < < 1.4 GeV. The data contain detailed in- 
formation on angular distributions to determine, at least 
qualitatively, the angular momentum content of the re- 
action amplitude. 

As is repeatedly mentioned in the paper, among differ- 
ent mechanisms which might contribute to the reaction, 
the excitation of the -D33 wave turns out to be most ev- 
ident in the data. According to our calculation, at the 
lab photon energy Ej = 1.1 GeV and above the coherent 
sum of £)33(1700) and 1)33(1940) can account for the ma- 
jority of the experimental results. We pay some attention 
to a qualitative discussion of the way in which the main 
features of the angular distributions determine some of 
the adjustable parameters of the model. 

At the same time, deviations from the single model 
predictions point to a nontrivial admixture of other reso- 
nance states. That these 'weaker' resonances should have 
positive parity is confirmed by visible forward-backward 
asymmetry in the cos 0,, distribution in the overall cm. 
frame (fifth column in Fig. [5]). Inclusion of ^33(1600), 
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P3i(1750), and ^35(1905), improves the quality of our 
fit. 

Furthermore, in order to maintain tire rather monotone 
rise demonstrated by the measured total cross section up 
to — 1.4 GeV we are forced to include the tails of the 
resonance states with higher masses. Taking into account 
the results obtained in ref. [sj, the resonance P33(1920) as 
well as the one-star resonance D33(1940) were included 
into the model. We allow their parameters to be varied, 
although the energy region covered by the measurements 
does not permit an unambiguous determination of their 
values. 

On the whole, our results, at least those corresponding 
to solution I, are rather close to the results obtained in 
ref. 0. We also see that a more rigorous study of the 
partial wave structure of tt"?] photoproduction requires 
detailed polarization data. Otherwise, the theory admits 
a variation of the resonance parameters over a relatively 
wide range. As an example, we give three versions of 
the fit which explain the features of the observables via 
different mechanisms. In particular, the characteristic 
shape of VF(O^) at higher energies, with a minimum at 
cos O^ = 0, may be described via different dynamical 
contents. In version III this shape is reproduced if we 
allow for a small admixture of the d-wave mode in the 77A 
decay of the 1)33(1700). In the other two cases the same 
effect was provided by the contribution of P33 and/or F35 
states. 

We also explored in considerable detail the photon 



asymmetry S, and found that its main features may well 
be interpreted in terms of the single D33 model. In this 
model the nontrivial value of E(7riV) is a simple conse- 
quence of "qp interaction in the final state. Inclusion of 
other resonances leads to rather small corrections. How- 
ever, as pointed out in Sect. lIII B| the behavior of Y,{'i]N) 
in the region of low rjN invariant energies is quite sen- 
sitive to the contribution of the P33 wave. In partic- 
ular, our solution II containing a large contribution of 
P33(1920) turns out to be inconsistent with the available 
data for j:{r]N). 

Further study of the reaction -> Tr^rjp may be aimed 
at those observables, in which the contributions of pos- 
itive parity resonances may be identified in a model in- 
dependent way, wherever possible. With new data for 
polarization observables one may expect considerable im- 
provement of accuracy in the determination of the res- 
onance parameters. In the past year some progress has 
been made towards that goal in refs. 15 1 and |16| . 
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